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(57) Abstract 

A method and system for determining 
marginal values (34) for perishable resources 
(31) expiring at a future time, for example, an 
airline seat, hotel room night, rental car day or 
the like, for use in a perishable resource revenue 
management system (MVS) (16). Departures, 
flight departiires and passenger value data are 
loaded into the MVS (16) from the airline 
revenue management system (50). These data 
and demand are constructed for storing the 
received information (51). Mar^ginal values (34) 
for each of the seats on the flight leg departures 
(31) are determined using an iterative function 
(52). Upon completion of the iterative function, 
tfie marginal vahies are stored back (53). 
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METHOD AND SYSTEM FOR DETERMINING MARGINAL VALUES 
FOR USE IN A REVENUE MANAGEMENT SYSTEM 

5 

Field of the Invention 

The present invention relates to a method and system for determining marginal values for 
use in a revenue management system and, in particular, to a method and system for determining 
marginal values for use in conjunction with an airline reservations system for controlling 
1 0 booking of seats on flight leg departxires. 
Background of the Invention 

Marginal values (or bid prices) are widely used in revenue management practice as an aid 
to guiding the setting of prices at which to sell or commit perishable resources. For instance, in 
the airline industry, marginal values can be used by an airline reservations system for granting 
1 5 or denying the sale of seats on one or more flight leg departures. 

Each marginal value represents the minimum price that should be accepted by a merchant 
to sell or commit an additional perishable resource based on a goal of maximizing profitability 
in the face of possibly uncertain demand. For a combination of perishable resources, the sum 
of individual marginal values represents the minimum acceptable price for the perishable 
20 resource combination. 

Perishable resources by definition are ones tiiat cannot be inventoried and share three 
common characteristics: perishability, "fixed" capacity and segmentability. Perishability means 
that each resource ages or becomes unavailable, and thus has no value, after a certain date, time 
or similar temporal event. "Fixed" capacity implies a high cost of adding an incremental unit 
25 such that capacity is regarded as static and unchanging. Segmentability refers to the ability to 
segment customers based on a willingness to pay using different rates and/or different purchase 
restrictions, such as the date of purchase relative to the date of use. Examples of perishable 
resources include airiine seats, hotel room nights, rental car days and similar products or services 
such as described in Weatherford & S.E. Bodily, A Taxonomy and Research Overview of 
30 Perishable-Asset Revenue Management: Yield Management, Overbooking, and Pricings 40 
Operations Research 5, pp. 831-44 (1992), the disclosure of which is incorporated herein by 
reference. 

One problem faced in the determination of marginal values is processing time. In the 
airline industr>', a single airline might fly possibly thousands of flight legs (a part of a flight 
35 consisting of a single take off and landing) each day. Each flight leg requires a marginal value 
which is itself dependent on the particular capacity of the airplane used and the known or 
expected demand for that flight leg as determined by the demand for all sequences of flight legs 
that carry passengers from their point of origin to their point of destination. The system 
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described herein can calculate marginal values for one departure day for a large airline in less 
than 1 5 minutes. However, prior art systems which account for these factors are constructed to 
calculate marginal values for all future departure days at one time, that is, during a single 
processing mn, and require substantially more than a day (and usually three days) (>f calculation 
time. 

Another problem faced is detennining whether the system has converged on a final 
solution. Prior art systems generally detennine marginal values using an itenitive process 
whereby candidate marginal values are refined during each successive iteration until a 
satisfactory set of values is obtained. One prior art approach places an arbitrary limit on the 
number of iterations performed. This limit is used to cut down on the processing time required 
and to avoid problems with the system oscillating between candidate solution sets. There is no 
assurance that the final "solution" has actually converged. Moreover, if the limit is set too high, 
the processing time required becomes unnecessarily long. Conversely, if the limit is set too low, 
the resulting set of marginal values might be far fi-om those that generate the maximum net 
revenue of the entire network. 

Yet another problem faced is how to treat dependencies between individual perishable 
resources. In the airline industry, a combination of flight legs (known as a flight path), such as 
occurs whenever a passenger books a series of coimecting flight legs, creates a dependency 
between each individual flight leg in the flight path, yet prior art systems generally overlook the 
effect that the flight path as a whole has on the individual marginal values. Experience with 
revenue management systems that take into account the passenger's flight path indiciite that there 
is substantial value in doing so, on the order of 1% to 2% of revenue compared with conventional 
leg-based systems. 

One further complication that has been recognized by the airline industr>', yet can be 
found in similar industries, is that of determining marginal values for flight leg departures in a 
multi-dimensional problem space. Prior art hotel industry approaches were merely concemed 
with a one-dimensional mapping of room nights based on each guest's length of stay. In the 
airline industry, though, marginal values must be determmed with a two-dimensional mapping 
of seats based on individual flight legs being associated with one or more fliglit paths and 
individual flight paths being associated with one or more flight legs. 

Therefore, there is a need for a method and system for determining optimal, system-wide 
net revenue maximizing marginal values for perishable resources, such as seats on flight leg 
departures, that does not require substantial processing time yet utilizes a convergence criteria 
for dynamically determining whether a substantially optimal set of marginal values has been 
obtained. Desirably, such a method and system must also factor in dependencies between 
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individual perishable resources, including the mapping of the determination in a multi- 
dimensional problem space. 

5 Summary of the Inventi n 

Accordingly, the present invention enables the above problems to substantially be 
overcome by providing a method and system for deteraiining marginal values for use in a 
revenue management system and a method and class library for adaptively determining client- 
specific marginal values in a generic problem space applicable to other industries, including 

10 hotels, rental car companies and other industries offering perishable resources. 

An embodiment of the present invention is a method and system using a marginal value 
system for determining marginal values for perishable resources expiring at a future time, for 
example, an airline seat, hotel room night, rental car day or the like. Data for the perishable 
resources and composite resources is loaded from a perishable resource revenue management 

1 5 system into the marginal value system. Internal data structures are constructed for linking each 
of the perishable resources to their associated composite resources and for linking each of the 
composite resources to their associated perishable resources. The marginal values for the 
perishable resources are determined using a continuous optimization function using 
interdependencies among the perishable resources and the composite resources in the internal 

20 data structures. The marginal values are stored from the marginal value system into the 
perishable resource revenue management system, 

A further embodiment of the present invention is a perishable resource revenue 
management system for granting and denying a sale of one or more perishable resources expiring 
at a futiue time depending on marginal values for each of such perishable resources received 

25 from a perishable resource revenue management system. The marginal values are determined 
by the steps of loading data for the perishable resources, composite resources and perishable 
resource values from a perishable resource revenue management system into a marginal value 
system; constmcting internal data stmctures for the perishable resources, composite resources 
and perishable resource values; determining marginal values for the perishable resources using 

30 a continuous optimization function; and storing the marginal values from the marginal value 
system into the perishable resource revenue management system. 

A further embodiment of the present invention is a method using a computer and a 
marginal value system for determining marginal values for perishable resources expiring at a 
future time. A locally optimal marginal value is evaluated for one of the perishable resources 

35 using a continuous optimization function dependent on the marginal values for the other 
perishable resources. The locally optimal marginal value is iteratively reevaluated xmtil a 
globally optimal marginal value is attained. 
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A further embodiment of the present invention is a method and a marginal value system 
for determining marginal values for seats on flight legs departing on a future date. Data for flight 
5 leg departures and flight path departures is loaded from an airline revenue management system 
into the marginal value system. Internal data structures are constructed for linking each of the 
flight leg departures to their associated flight path departures and for Hnking each of the flight 
path departures to their associated flight leg departures. The marginal values for seats on the 
flight leg departures are determined using a continuous optimization function using 
1 0 interdependencies among the flight leg departures and the flight path departures in the internal 
data structures. The marginal values are stored fix)m the marginal value system into the airline 
revenue management system. 

A further embodiment of the present invention is a method using a computer and a 
marginal value system for determining marginal values for seats on flight legs departing on a 
15 future date. A locally optimal marginal value is evaluated for a seat on one of the flight leg 
departures using a continuous optimization function dependent on the marginal vdues for a seat 
on the other flight leg departures. The locally optimal marginal value is iteratively reevaluated 
until a globally optimal marginal value is attained. 

Still other embodiments of the present invention will become readily apparent to those 
20 skilled in the art from the following detailed description, wherein is shown and described only 
embodiments of the invention by way of illustration of the best modes contemplated for carrying 
out the invention. As will be realized, the invention is capable of other and different 
embodiments and its several details are capable of modification in various obvious respects, all 
without departing from the spirit and scope of the present invention. Accordingly, the drawings 
25 and detailed description are to be regarded as illustrative in nature and not as restrictive. 
Brief Description of the Drawings 

FIG. 1 is, by way of example, an outline map of the topology for an airline network; 
FIG. 2 A is a functional block diagram of an airline reservations booking system; 
FIG. 2B is a flow diagram of a method for using marginal values in the airline 
30 reservations booking system of FIG. 2 A; 

FIG. 3 shows data structures used by the marginal value system of FIG. 2A ; 
FIG. 4 is a flow diagram of a method for determining marginal values for use in the airline 
reservations system of FIG. 2 A; 

FIG. 5 is a flow diagram of an iterative function for determining marginal values for flight 
35 leg departures; 

FIG. 6 is a flow diagram of a supply-demand balance optimization function used by the 
function of FIG. 5; 
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FIGS. 7A and 7B shoAv graphs of flight leg departure demand curves generated by the 
flight leg departure demand determination function of FIG. 8; 
5 FIG. 8 is a flow diagram of a flight leg departure demand determination function; 

FIG. 9 is a flow diagram of a flight path departure demand determination function; 

FIG. 1 0 is a flow diagram of a flight path departure demand search function; 

FIG. 1 1 A shows a graph of a flight path departure demand curve used by the search 
function of FIG. 10; 

10 FIG. 1 IB is a table showing, by way of example, the results of a demand curve 

construction; 

FIG. 12 is a flow diagram of an expected marginal seat revenue (EMSR) optimization 
function used by the function of FIG. 5; 

FIG. 13 shows a graph of a flight leg departure protection level curve generated by the 
1 5 EMSR optimization function of FIG. 1 2; 

FIG. 1 4A is a flow diagram of a demand point list creation function; 

FIG. 148 is a data structure used by the EMSR optimization function of FIG. 12; 

FIG. 14C is a table showing, by way of example, the calculation of protection levels; 

FIG. 1 5 is a flow diagram of a protection level determination function; 
20 FIG. 1 6 is a class diagram of a class library for a method for determining marginal values 

for use in a revenue management system; and 

FIG. 17 is a class diagram of an instantiation of the class library of FIG. 16 for an airline 
reservations system; 

FIG. 18 is a block diagram of a method for determining generic marginal values for use 
25 in a client-specific revenue management system; 

FIG. 1 9 is a flow diagram of an iterative function for determining generic marginal values 
for perishable resources; 

FIG. 20 is a flow diagram of a supply-demand optimization function used by the function 
of no. 19; 

30 FIG. 2 1 is a flow diagram of a resource demand determination function; 

FIG. 22 is a flow diagram of a composite resource demand determination function; 
FIG. 23 is a flow diagram of a composite resource demand search function; 
FIG. 24 is a flow diagram of a non-deterministic optimization function used by the 
function of FIG. 19; 

35 FIG. 25 is a flow diagram of a demand point list creation function; and 

FIG. 26 is a flow diagram of a protection level determination function. 
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Detailed Descripti n 
I. Glossary of Terms 

5 "Flight" refers to a sequence of regularly scheduled aircraft take oflFs and landings 

designated by a flight number. A flight begins at an origin airport, can make one or more stops 
and terminates at a destination airport. "Flight leg" refers to one part of a flight consisting of 
a single take off and landing. 

"Flight leg departure" refers to a flight leg flown on a specific date (called the departure 

10 date). 

"Flight Path" or "itinerary" refers to one or more connecting flight legs flown by an 
individual passenger from an origin airport to a destination airport. 

"Flight path departure" refers to a flight path flown on the departure date. 

"Airline network" refers to the collection of flight legs flown by an airline. 
15 "Booking class" refers to a code used in reserving or booking a seat on a flight leg 

departure and is related to a ticket's price and conditions of purchase. 

"Demand forecast" represents a prediction of the number of seats that passengers will 
request to book for a future flight path departure specified by statistical parameters, such as mean 
and variance. It also can include forecasts of cancellations and no-shows. 
20 "Passenger value" represents what an airline receives for a passenger flying on a given 

flight path minus any associated direct costs, such as taxes, commissions and carrying costs, for 
example, the meal served and the additional fuel needed to carry that passenger. 

"Resource" or "perishable resource" refers to a single unit of product or service that a 
merchant can sell or commit for sale. A flight leg departure is an example of a resource. 
25 "Composite resource" refers to a collection of resources that a customer typically will 

purchase as a bundle at a single price. For example, a flight path departure is a composite 
resource made up of flight leg departures which are the resources. 

IL Overview 

30 

A. Airline Network Topology 

FIG. 1 shows, by way of example, an outline map of the topology for a simple airline 
network 9. The auiine network 9 comprises a plurality of airports as exemplified below in Table 
I. Each of the airports is intercoimected to at least one other airport by one or more flight legs 
35 10. The airline network 9 is configured with a single hub located at JFK. A representative flight 
number of each flight leg 1 0 is labeled on the line connecting each pair of airports. For instance, 
the flight leg 1 0 originating at LAX and terminating at JFK is flight 5 1 . 
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Table I 

C de Description 

LAX Los Angeles International Airport (Los Angeles, CA) 

ORD O'Hare International Airport (Chicago, IL) 

BOS Logan International Airport (Boston, MA) 

JFK John F. Kennedy International Airport (New York, NY) 

PHL Philadelphia International Airport (Philadelphia, PA) 

LHR London Heathrow Airport (London, U.K.) 

CDG Charles de Gaulle Airport (Paris, France) 



B. Airiine Reservations Booking System 

FIG. 2A is a functional block diagram of an airline reservations booking system 1 1 , used, 
by way of example, for granting or denying the booking of seats on flight leg departures based 
on marginal values. The booking system 1 1 comprises five components: a travel agent terminal 
12, a global distribution system (GDS) 13, an airline reservations system 14, an airline revenue 
management system 1 5 and a marginal value system (MVS) 16. Although only one of each type 
of component is shown, a typical airline reservations booking system 1 1 would comprise 
multiple travel agent terminals 12, GDSs 13, airline reservations systems 14, airline revenue 
management systems 1 5 and marginal value systems 16. 

Each GDS 13, airline reservations system 14, airline revenue management system 15 and 
marginal value system 16 can be embodied as a general purpose progmmmed computer system 
(not shown) that operates on a set of instructions in the form of a computer program. Such a 
general purpose programmed computer system contains a central processing unit, volatile 
memory, secondary storage and other components conventional in the art, including a cathode 
ray tube (CRT) or other type of display and a keyboard entry device. 

Each travel agent terminal 12 is a typical computer terminal, including a cathode ray tube 
(CRT) or other type of display and a keyboard entry device (not shown) and is used by travel 
agents for remotely accessing the GDS 1 3 over dedicated data lines. In FIG. 2 A, only one travel 
agent terminal 12 is shown, although there can be a multiplicity of these terminals connected to 
the GDS 13. 

The travel agent terminal 12 is used for submitting availability and seat booking requests 
to the GDS 1 3 and for receiving data in reply from the GDS 1 3. An availability request queries 
whether a seat is available to take a passenger &om an origin airport to a final destination airport 
over possibly different flight paths flown by possibly multiple airlines. A seat booking request 
specifically asks for a seat reservation. There are multiple airline reservations systems 14 
associated with the GDS 13. To determine the aveiilability of seats on the flight legs 10 flown 
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by each airline, the GDS 1 3 sends availability and booking requests to each airline reservations 
system 14 involved in the travel agent's request. Information about the flight leg departures 10 
in the airline network 9 (shown in FIG. 1) is managed and stored in the airline reservations 
system 14. 

A response to an availability request is presented as a series of coded booking classes 
paired with an availability indicator. As is conventional in the art, a series of eight booking 
classes are used, although the exact nimiber and combination is arbitraiy and could vaiy 
according to the application. An example of a response is "F5C8 Y9B9M5H4Q2V0," where each 
alpha and numeric character pair indicates a class of service and the approximate number of seats 
available. For instance, " Y9" indicates that there are 9 or more seats available on the particular 
flight path departure for "Y" class, "M5" indicates that there are 5 seats available for "M" class 
and **V0" indicates that no seats are available for "V" class. The travel agent uses a booking 
class to submit a booking request to the airline reservations system 14 through the GDS 13. 

The GDSs 13 and the airline reservations systems 14 are shown as separate blocks in FIG. 
2A, although in practice these systems are sometimes operated on the same computer. The 
airline reservations system 14 can also perform many of the functions of the GDS 13, although 
it focuses on the services of a single airline. Participation by an airline in the GDS 1 3 is optional, 
although most airlines participate in most of the major GDSs in the world. In the United States, 
examples of GDSs 13 are SABRE, APOLLO, WORLDSPAN and SYSTEM ONE, In Europe, 
two consortium GDSs 13 include AMADEUS and GALILEO. 

The airline reservations system 14 accepts or denies the booking requests that it receives 
from the GDS 13 or a locally-connected terminal (not shown) based on various factors, including 
mai^inal values received from the airiine revenue management system 1 5. Generally, the airline 
revenue management system 15 is also owned by an individual airline and operates on a separate 
computer from the airline reservations system 14. Its purpose is to calculate the seat availability 
according to booking class and corresponding to customer origin and destination jtnd departure 
date and to process booking requests using stored marginal values. 

Part of the airline revenue management system 15 is the marginal value sj'slem (MVS) 

16. In the described embodiment, the MVS 16 is a computer program written in the C+h- 

pogramming language that is executed by the airline revenue management system 1 5, although 

other programming languages are feasible. It is also envisioned that the MVS 1 6 can be a special 

purpose computer formed on an integrated circuit chip (or set of chips) or executable computer 

code burned into a read-only memory (ROM) chip that can be read in by conventional means or 

as microcode. An exemplary example of a marginal value system 16 is the marginal value 

TM 

system product licensed under the trademark of Marginal Value Engine manufactured by 
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Decision Focus Incorporated, Mountain View, CA. Marginal Value Engine and MVE are 
trademarks owned by Decision Focus Incorporated. 

5 The purpose of the MVS 1 6 is to determine system- wide optimal marginal values for use 

by the airline revenue management system 15. The MVS 16 periodically receives a demand 
forecast, passenger value and the supply of seats remaining to be booked for selected future flight 
leg departures as inputs from the airline revenue management system 15. This same data could 
also be retrieved from a storage device (not shown) commonly accessible by both the airline 

1 0 revenue management system 15 and the marginal value system 16. A demand forecast predicts 
the number of bookings that will be or is expected to be made between now and the departure 
date expressed in terms of statistical parameters, such as mean and variance. The passenger 
value is a function of the booking class, point of sale and the flight path based on origin and 
destination. The supply of seats remaining to be booked is calculated from the capacity of the 

1 5 airplane assigned to the flight leg departure minus the seats already booked plus an overbooking 
factor to account for cancellations and no-shows. The number of flight leg departures that can 
be solved at one time is limited only by the computer memory available. A very large airline 
network 9 can be solved on a large computer. 

In response, the MVS 16 sends substantially optimal marginal values for the flight leg 

20 departures selected to the airline revenue management system 1 5 for immediate use or storage 
or stores them directly into the commonly accessible storage device (not shown). 

In the described embodiment, the MVS 16 is a computer program vmtten in the C-h- 
object-oriented programming language. It is preferably structured as a class library as further 
described hereinbelow. However, any suitable programming language can be employed and the 

25 MVS 1 6 can be embodied on any suitable computer system configuration. 

C Method for Using Marginal Values 

FIG. 2B is a flow diagram of a method for using marginal values in the airline 
reservations system 14 of FIG. 2 A. Marginal values determined by the MVS 1 6 are used by the 

30 airline reservations system 14 in granting or denying the sale of seats on one or more of the flight 
legs 10 departing on some future date. It is feasible that the marginal values determined by the 
MVS 16 could be used in different ways than described herein by airline reservations systems 
1 4 and airline revenue managemait systems 1 5. Moreover, it is feasible that the event triggering 
the use of marginal values is a time, week, month, year or other future temporal event instead of 

35 or in addition to a date. 

A booking request is received by the airline reservations system 14 for a flight path 
departure from the travel agent terminal 12 through the GDS 13 (block 18). The airline 
reservations system 14 determines a summation of the marginal values for the flight leg 
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departures contained in the flight path departure (block 1 9) to establish a total marginal value 
for the origin-to-destination flight path departure. The airline reservations system also 

5 determines the net revenue for the same flight path departure (block 20). The summation of the 
marginal values is compared to the net revenue and if the net revenue exceeds the summation of 
the marginal values for the fli^t path departure (block 2 1 ), an acceptance of the booking request 
is sent to the travel agent terminal 12 through the GDS 13 (block 23). Otherwise, if the net 
revenue is less than the summation of the marginal values (block 2 1 ), it is not profitable to accept 

10 the booking request and a denial of the booking request is sent to the travel agent terminal 12 
through the GDS 13 (block 22), 

IIL Airline-Specific Marginal Value Determination Method 
A. Data Structures 

1 5 FIG. 3 shows the data structures used by the marginal value system 1 6. Three major data 

structures are used: a flight leg departures list 30. a flight path departures list 35 ;md a demand 
records list 38. 

The data structure for the flight leg departures list 30 comprises flight leg departing 
records 31 containing fields for storing a departure date 43, flight number 32, origin and 

20 destination airports 33, marginal value 34 and supply of seats 44. The data stricture for the 
flight path departures list 35 comprises flight path departure records 36 containing fields for 
storing a departure date 46 and a flight path definition 37. The data structure for the demand 
records list 38 comprises demand records 39 containing fields for storing a booking class price 
40, mean 41 and variance 42. In the described embodiment, a departure date 43 is specified, 

25 although a departure time or other future temporal event is feasible. 

Each flight in the airline network 9 is divided into its respective flight legs and the 
corresponding flight leg departures 10 are stored in the flight leg departures list 30. The flight 
number 32 field links each flight leg departure record 31 to its respective flight in the airline 
network 9, The origin and destination 33 field indicates the part of the flight to which the flight 

30 leg 1 0 corresponds. The departure date 43, marginal value 34 and supply of seats 44 fields store 
the departure date, marginal value and seat supply for the flight leg departure 10. For 
convenience, the flight leg departure records 31 will be referenced whenever a flight leg 
departure is refen^d to hereinbelow. 

The data stmcture for the flight path departures list 35 stores the actual or forecasted flight 

35 path departures through the airline network 9. The departure date field 46 stores die departure 
date of the first flight leg 1 0 contained in the flight path departure. The flight path definition 37 
field defines the flight path departure in terms of origin, intermediate and destiruition airports 
with the connecting flight number indicated between each of the airports. For convenience, the 
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flight path departure records 36 will be referenced whenever a flight path departure is referred 
to hereinbelow. 

5 A set of pointers 45 interconnect the flight leg depaitiues list 30 and the flight path 

departures list 35 to form a two-dimensional mapping of interdependencies between the two data 
structuies. Each flight leg departure record 31 has a list of pointers 45 pointing to all flight path 
departure records 36 containing the flight leg departure 10 stored in the flight path deflnition 
field 37 of the flight path departure record 36. Similarly, each flight path departure record 36 

10 has a list of pointers 45 pointing to all flight leg departure records 3 1 containing the flight leg 
departure 10 stored in the flight path definition field 37 of each flight path departure record 36. 
For simplicity, the data structures for these pointer lists arc not shown in FIG. 3. 

TTie data structure for the demand records list 38 stores the different fares and statistical 
parameters from the demand forecast for the corresponding flight path departure. The fares are 

1 5 stored in the booking class price 40 field in decreasing order. Two statistical parameters, mean 
and variance, are stored in the mean 41 and variance 42 fields, although further statistical 
parameters are feasible. Each flight path record 36 points to one demand records list 38 which 
contains as many demand records 39 for as many fares as required. For convenience, the 
booking class price 40 field will be referenced whenever a particular booking class price is 

20 referred to hereinbelow, 

B. Method for Determining Marginal Values 

In the described embodiment, different optimization functions can be employed in 
determining marginal values for seats on the flight leg departures based on assumptions that the 

25 demand is either deterministic or non-deterministic. It is critical that the optimization function 
employed generate a continuous demand curve, preferably as a function of marginal value 
(expressed in dollars) and supply (expressed in seats on the flight leg departure). Two 
optimization functions, a supply-demand balance and an expected marginal seat revenue 
(EMSR), arc described herein; however, other optimization functions are feasible depending 

30 upon the characteristics of the particular embodiment, such as booking patterns, marginal value 
updating frequency and events that trigger reoptimization. 

The basic control structime used to iteratively determine marginal values is identical 
regardless of the optimization function employed. This generic control structure is shown in 
FIGS. 4 and 5. A supply-demand balance optimization fimction is depicted in FIGS. 6 through 

35 1 1 B. An EMSR optimization function is depicted in FIGS. 1 2 through 1 5. The iterative section 
common to both optimization functions will now be described. 

FIG. 4 is a flow diagram of a method using the MVS 1 6 for determining marginal values 
for use in the airline reservations system 14 of FIG. 2A. Information regarding flight leg 
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departures, flight path departures and passenger value data are loaded into the M\^S 16 from the 
airline revenue management system 15 (block 50). The internal data structures shown in FIG. 

5 3, including the flight leg departures list 30, the flight path departures list 35 and the demand 
records list 38, are constructed for storing the received information (block 5 1 ). Marginal values 
34 for each of the seats on the flight leg departures 3 1 are determined using an iterative function 
(block 52) as further described hereinbelow in FIG. 5. Upon completion of the itenitive function, 
the marginal values are stored back to the airline revenue management system 15 (block 53). 

10 FIG. 5 is a flow diagram of an iterative function (block 52) for determining marginal 

values 34 for the flight leg departures 31 as stored (block 51) in the data structures of FIG. 3. 
This function is structured as a pair of nested control loops that sequentially determine a 
marginal value 34 for each flight leg departure 3 1 using a Gauss-Seidel iterative approach. The 
fligjit leg departures 3 1 and the flight path departures 36 are interdependent upon each other for 

15 purposes of marginal value determination. The outer control loop (blocks 61-67) repeats the 
determination of marginal values 34 for all of the flight leg departures 31 until a convergence 
criterion is satisfied. In the described embodiment, convergence is generally attained in as few 
as 3-4 iterations, but almost always in fewer than 1 0 iterations resulting in a substantial decrease 
in processing time over the prior art. The inner control loop (blocks 62-66) sequentially 

20 determines a marginal value 34 for each of the flight leg departures 3 1 . 

The marginal value 34 for each of the flight leg departures 3 1 is initialized to 0 (block 60). 
A local variable maxjliffxs set to 0 (block 61). This variable stores the largest difference 
between current versus previous marginal value 34 out of all marginal values 34 determined 
during the current iteration of the outer control loop (blocks 61-67) and is used to determine 

25 convergence. A locally optimal marginal value 34 is evaluated for a seat on one of the flight leg 
departures 31 using a continuous optimization function (block 62). In the described 
embodiment, two continuous optimization functions, supply-demand balance and EMSR, are 
used as further described hereinbelow in FIGS. 6-1 IB and FIGS. 12-15, respectively. 

The continuous optimization fimction returns a marginal value 34 for a seat on the current 

30 flight leg departure 31. The magnitude of change in the locally optimal marginal value 34 
(obtained in block 62) is determined (block 63) as expressed by the equation: 

diff = - MVi_^\ 

35 where diff \s a local variable and MV^ and A^/.y are the current and previous marginal values 
34, respectively, for the current flight leg departure 3 1 . D/^measures the change which occiured 
between the current marginal value 34 as calculated during the present iteration of the outer 
control loop (blocks 61-67) versus the previous marginal values 34 as calculated during the 
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previous iteration, \fdiff \s greater than max c//j^(block 64), max diff 'is set to diff (block 65) 
since rff^reflects the biggest change in marginal value that has occurred up through this iteration 

5 of the inner control loop (blocks 62-66). Otherwise, mar ^i/f remains unchanged (block 64). 

If other flight leg departures 31 remain (block 66), that is, must be evaluated for a locally 
optimal marginal value 34, control returns to the top of the inner control loop (block 62-66) for 
farther processing. Otherwise, if the last flight leg departure has been evaluated (block 66), the 
inner control loop (block 62-66) is exited. 

10 In the described embodiment, a convergence criterion is used in determining whether 

globally optimal marginal values 34 have been found. Consequently, the outer control loop 
(blocks 61-67) is performed for an indefinite number of iterations until the maximum difference 
(as stored in the local variable max diff) between the current and previous marginal values 34 
out of all of the flight leg departures 3 1 is less than the convergence criterion. In the described 

15 embodiment, a convergence criterion of 0.02 (corresponding to 2 cents) is used, although any 
other desired convergence criterion value could be used. 

Thus, if maxji\ff \% less than convergence (block 67), control returns to the top of the 
outer control loop (block 61-67) for farther processing. Otherwise, a globally optimal marginal 
value 34 has been obtained for each of the flight leg departures 31 (block 67) and the fanction 

20 returns. 

A maximum iteration linut is used to prevent the fanction from rurming endlessly in the 
case of an error. If the iteration count exceeds the maximum iteration limit (block 68), the 
system will report the error condition of non-convergence (block 69) and return. In the described 
embodiment, a maximum iteration limit of 1 00 is used, although any other desired iteration limit 
25 could be used. 



C. Supply-Demand Balance Optimization Function 

FIG. 6 is a flow diagram of a supply-demand balance optimization function (block 62) as 
used by the iterative fanction of FIG. 5. Its purpose is to find a locally optimal marginal value 
30 34 based on the assumption that the demand for seats is deterministic and that the optimal 
marginal value is that which substantially balances supply and demand. 

The marginal value 34 for the current flight leg departure 3 1 is initialized to 0 (block 70). 
The demand for the current flight leg departure 3 1 is determined using a demand determination 
function that creates a demand curve as a function of marginal value 34 (shown in FIGS. 7A and 
• 35 7B) as farther described hereinbelow in FIG. 8, The difference between the supply of seats 44 
available on the current flight leg departure 31 and the corresponding demand for seats (as 
returned by the determination function of FIG. 8) is determined (block 71). If the difference of 
supply 44 minus demand is greater than 0 (block 72), there are more seats available on the 
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current flight leg departure 31 than in demand, that is, supply 44 exceeds demand, and the 
marginal value 34 remains at 0. The function then returns. Otherwise, if the difference of supply 
44 minus demand is less than or equal to 0 (block 72), the demand determination function of 
FIG. 8 is repeatedly called and the difference between supply 44 and demand taken until a 
marginal value 34 that makes supply 44 and demand equal is found (block 73). The function 
then returns. The operation of this function will now be described with reference to FIGS. 7A 
and 7B which show graphs of flight leg departure demand curves 8 1 and 85, respectively. 

The graphs of FIGS. 7A and 7B illustrate the two cases selected by the decision block 72 
of FIG. 6. The x-axis represents marginal values expressed in dollars and the y-axis represents 
seats on the current flight leg departure 3 1 . The supply of seats 44 (as indicated by line 80 in 
FIG. 7A and by line 84 in FIG. 7B) is provided by the airline revenue management system 15 
as an input parameter. The demand curve (as indicated by line 81 in FIG. 7A and by line 85 in 
FIG. 73) is generated by the flight leg departure demand determination function <:>f FIG. 8. 

By way of example, assume that a marginal value 34 for the ORD-JFK flight leg departure 
3 1 is being determined with a supply of seats 44 for this flight leg departure equalling 1 0. In 
FIG. 7A, the situation where the difference between supply 44 minus demand is greater than 0 
is shown. At point 79, the marginal value 34 equals 0, supply 44 equals 10 and demand equals 
8. Thus, the difference 82 between supply 80 and demand 81 equals 2 and is therefore greater 
than 0. This situation indicates that the supply 44 exceeds the demand and any non-zero 
marginal value is therefore meaningless. In such cases, a marginal value 34 equaling 0 will be 
retimied. 

In FIG 7B, the situation where the difference between supply 44 minus demand is less 
than or equal to 0 is shown. At point 83, the marginal value 34 equals 0, supply 44 (equals 1 0 and 
demand equals 1 4. Thus, the difference 86 between supply 84 and demand 85 equals -4 and is 
therefore less than 0. This situation indicates that the demand exceeds the supply 44 and a non- 
zero marginal value is required. Thus, the intersection 87 between supply 84 and demand 85 
must be found. In the described embodiment, a Van Wijngaarden-Dekker-Brent search method, 
such as described in W. H. Press et al,. Numerical Recipes in C-The Art of Scientific Computings 
pp. 267-69, Cambridge University Press (1988), the disclosure of which is incorporated herein 
by reference, is employed to search the demand curve 85 (block 73) using the flight leg departure 
demand detemiination function of FIG. 8 as an input parameter along with the supply 44. This 
function is used to generate candidate demand points 88A, 88B along the demand curve 85 as 
the Van Wijngaarden-Dekker-Brent search method searches for the intersection 87. In effect, 
the search moves along the x-axis using a candidate marginal value 34 at each candidate demand 
point 88A, 88B and the demand for seats at each candidate demand point 88A, 88B is compared 
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to the supply 44. Once the intersection 87 is found, the marginal value 34, here, equalling $386, 
is obtained as a point 89 along the x-axis. 

5 FIG. 8 is a flow diagram of a flight leg departure demand determination function (blocks 

71 and 73) as used by the supply-demand optimization function of FIG. 6. Its purpose is to 
generate the demand curve 81, 85 for the current flight leg departure 31 as a function of a 
candidate marginal value 34. The demand for the current flight leg departure 31 is initialized to 
0 (block 90). The demand for a flight path departtire 35 containing the current flight leg 

1 0 departure 3 1 is determined by a flight path departure demand determination function (block 9 1 ) 
as further described hereinbelow in FIG. 9. The value this function returns represents the 
demand for the flight path departure 36 at the candidate marginal value 34 for the current flight 
leg departure 31 and the current marginal values 34 for all other flight leg departures 31 
contained in the flight path departure 36. The demand for the flight path departure 36 is added 

15 to the demand for the current flight leg departure 31 (block 92). If there are more flight path 
departures 36 that contain the current flight leg departure 3 1 (block 93), the control loop (blocks 
91-93) is repeated until all flight path departures 36 that contain the current flight leg departure 
31 have been processed. In effect, the control loop (blocks 91-93) follows every pointer 45 
pointing from a flight leg departure record 31 to a flight path departure record 36 in a one-to- 

20 many relationship. When no more flight path departures 36 remain (block 93), the function 
returns. 

FIG. 9 is a flow diagram of a flight path departure demand determination function (block 
91) as used by the function of FIG. 8. Its purpose is to determine the demand for the current 
flight path departure 36 using the total of the marginal values 34 of the flight leg departures 3 1 

25 contained in the flight path departure 36. The total marginal value for the current flight path 
departure 36 is initialized to 0 (block 100). The marginal value 34 associated with one of the 
flight leg departures 31 contained in the current flight path departure 36 is added to the total 
marginal value (block 101). If other flight leg departures 31 contained in the flight path 
departure 36 remain (block 1 02), control returns to the top of the control loop (blocks 101-1 02). 

30 In effect, the control loop (blocks 1 01 - 1 02) follows every pointer 45 pointing from a flight path 
departure record 36 to a flight leg departure record 31 in a one-to-many relationship. When the 
marginal values 34 for all of the flight leg departures 3 1 contained in the flight path departure 
36 have been added to the total marginal value (block 102), the demand for the flight path 
dqjarture 36 is determined by a flight path departure demand search function (block 103) using 

35 the total marginal value as further described hereinbelow in FIG. 10. The function then returns. 

FIG. 10 is a flow diagram of the flight path departure search function (block 103) used 
by the fiinction of FIG. 9. Its purpose is to find points on a flight path departure demand curve 
121 (shown in FIG. 1 1 A) that bracket the supply and to perfonm a linear interpolation of the 
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demand curve 121 for the flight path departure 36 as a function of the total marginal value. The 
search and linear interpolation will be explained with reference to FIG. 1 1 A which shows a 
graph of the flight path departure demand curve 121. This curve is represented by an additional 
data structure (not shown) used by the supply-demand balance optimization function (shown in 
FIG- 6) and is created during the construction of the internal data structures in block 5 1 of FIG. 
4. 

In the graph of FIG, II A, the x-axis represents the price of the current flight path 
departure 36 expressed in dollars and the y-axis represents the demand expressed in seats. Each 
flight path departure record 36 points to an associated flight path departure demand curve 121 
which is constructed fix)ra the demand record Ust 38. 

The construction, by way of example, of the flight path departure demand curve 121 is 
accomplished by the follov^ng procedure. It is critically important that the demimd curve 121 
be continuous as guaranteed by the following procedure. Starting with the first demand record 
39 on the demand record list 38, that is, with the demand record 39 storing the highest booking 
class price 40, points of the form (price, demand) are added to the demand curve list 38. The 
first point in the example is: 

(5CPi, 0) 

where BCP| is the booking class price 40 of the first demand record 39 
the demand curve 12 1 of the form: 

(BCP^y cumulative jnean) 

where BCP^ is the booking class price 40 of the successive demand records / and 
cumulative jnean is the accumulation of the mean demand 41 using the foraiula: 

cumulative jnean += mean^^ 

where mean^^i is the mean demand of the demand record 39 before /. After reading the last 
demand record 39, one more point is added: 

(0, cumulative jnean) 

where cumulative mean includes the mean demand of the last demand record 39 in the demand 
record list 38. 



Points are inserted onto 



(2) 
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Once the demand curve 121 is constructed from points in the demand record list 140, 
additional points are inserted between each of the existing points in the demand curve 121, The 
form of the inserted point is: 

(price^y demandj) 



where pricey and demand^ are given by the fomiulas: 

price^ = price^^^ + (price^^^ - price^^^) x IF 

and 

demand^ - demand j + {demand^^^ - demand^^^ x IF 

where price and demand are the price and demand of the point on the demand curve 
following the inserted point and price^_j and demand-_j are the price and demand of the point 
on the demand curve preceding the inserted point. IF is a parameter called the interpolation 
factor. In the described embodiment, the value of the interpolation factor is 0.9, but any other 
value greater than 0.5 and less than 1 .0 could be used depending on the needs of the application. 

FIG. 11 B is a table showing, by way of example, the results of the demand curve 
construction using the demand records 39 (shown in FIG. 3). The odd numbered points (120A, 
120C, 120E, 120F, 120G, 1201, 120K) were inserted using formula (2). For example, point 120E 
has a price of $1 50 which comes from the third demand record 39 in the demand record lisl 38 
(having a booking class price 40 of $150) and a demand of 2 + 5 = 7 which is the accumulated 
demand of the first two demand records 39 in the demand record list 38. TTie even numbered 
points (120B, 120D, 120F, 120H, 120J) were calculated with formulas (4), (5) and (6). For 
example, point 120D is calculated as follows; 

pricei20D = 150 ^ ((400 - 150) x 0.9) = 375.0 

and 

demandi2QD = 2 + ((7 - 2) x 0.9) = 6.5 
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The points 120A, 120B, 120C, 120D, 120E, 120F, 120G, I20H, 1201, 120J, 120K in FIG. IIB 
are graphed as shown in FIG. 1 1 A. The prices are shown on the x-axis and are labelled 1 19A, 
1 19B, 1 190, 1 19D, 1 19E, 1 19F, 1 19G, 1 19H, 1 191, 1 19J, 1 19K. 

The graph of FIG, 1 1 A also illustrates the search and linear interpolation of the flight path 
departure demand curve 121 performed by the flight path departure search function of FIG. 10. 
The formula for interpolated demand is: 

(7W - pricey) 

interpolatedjiemand = demand — 

' (pricci ^ - pricey) 

^ {pricBi ^-TMV) 

+ demand : 

(price^ ^ - pricey) 

where price ^ and demand^ are the price and demand of the point with a price greater than the total 
marginal value, 7MF, and price ^-^j and demand^^j are the price and demand of the point with a 
price less than or equal to TMV. 

The fimction depicted in FIG. 1 0 starts with a search starting from 0 for the first demand 
point with a price greater than the total marginal value (blocks 1 10). If no such demand point 
is found (block 1 1 1), the function terminates and returns a value of 0 (1 13). FIG. 1 1 A shows, 
by way of example, a total marginal value equalling $260, as indicated at point 122. The first 
point with a price greater than $260 is (375, 6.5). The point (375, 6.5) and the previous point 
(1 50, 7.0) are interpolated (block 112) by applying formula (5) as follows: 

6,5 K mJLm . 7.0 X (375 - 260) ^ ^ 
(375 - 150) (375 - 150) 

The interpolation is shown graphically as the intersection point 124 of the line 123 ;ind the total 
marginal value 122 between the points 120d and 120e. The resulting demand value 125 is 
calculated and the function returns. 

D. Expected Marginal Seat Revenue (EMSR) Optimization Function 

FIG. 12 is a flow diagram of an EMSR optimization function (block 62) as alternatively 
used by the iterative function of FIG. 5. Its purpose is to find a locally optimal marginal value 
34 based on the assumption that the demand for seats is uncertain, that is, non-deterministic. 
This assumption is commonly used in leg-based airline yield management systems, such as 
described in P. Belobaba, Application of a Probabilistic Decision Model to Airline Seat 
Management Control, 37 Operations Research, pp. 1 83-97 ( 1 989) and RD. WoUmer, An Airline 
Seat Management Model for a Single Leg Route When Lower Fare Classes Book First, 40 
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Operations Research, pp. 26-37 (1992), the subject matters of which are incorporated herein by 
reference. 

The EMSR approach uses values known as protection levels in a similar way that demand 
is used in the supply-demand balance function. The EMSR method also uses a new data 
stmcture called a demand point list 140 (shown in FIG. I4B) and the concept of an accumulated 
demand point which will be described in the discussion of FIG. 12. The EMSR optimization 
function (shown in FIG. 12) will now be described. 

A demand point list 140 (shown in FIG. 14B) is created for the current flight leg departure 
3 1 using a function (block 130) further described hereinbelow in FIG, 14A. 

An accumulated demand point is initialized in block 1 3 1. An accumulated demand point 
consists of six local variables: accumulated mean 152, accumulated net leg revenue 153, 
accumulated variance 154, average net leg revenue 155, protection level 156 and next net leg 
revenue 157, several values of which are shown in FIG. 14C. The accumulated demand point 
is initialized by setting the accumulated mean 152, accumulated net leg revenue 153 and 
accumulated variance 154 to 0 (block 131). These variables are used for storing accumulated 
totals derived from the demand points in the demand point list 140 for the current flight leg 
departure 30. The average net leg revenue 1 55 and the protection level 156 are also initialized 
to 0. The initial value of the next net leg revenue 157 is set to the value of the net leg revenue 
142 of the first demand point 141 A in the demand point list 140. 

The first accumulated demand point 161 A is plotted on the graph shown in FIG. 13 using 
the next net leg revenue 157 on the x-axis and the protection level 156 on the y-axis. The 
remaining accumulated demand points I61B, 161C, 161D are iteratively determined using a 
control loop (blocks 132-137) that sequentially processes each demand point in the demand point 
list 140 starling with the point with the highest average net leg revenue 142. During each 
successive iteration of the control loop (blocks 132-137), the accumulated mean 152, 
accumulated net leg revenue 153 and accumulated variance are updated (block 132) based on 
the next demand point in the demand point list 140. The updated accumulated mean 152 is 
determined as expressed by the equation: 



where mean^ is the mean 143 for the next demand point in the demand point list 140. The 
updated accumulated net leg revenue (NLR) 1 53 is determined as expressed by the equation: 



accumulated mean mean: 



(8) 



accumulated NLR += NLR: x mean: 



(9) 
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10 



where NLR^ is the net leg revenue 142 and mean^ is the mean 143 for the next demand point in 
the demand point list 140. The updated accumulated variance is determined as expressed by the 
equation: 

accumulated_yariance += var^ ^2q^ 

where var^ is the variance 1 44 for the next demand point in the demand point list 1 40. 

The average net leg revenue (NLR) 155 is determined by dividing the accimiulated net 
leg revenue 153 by the accumulated mean 152 (block 132) as expressed by the equation: 

xrr D accumulated NLR /i i \ 

average _NLR = = 

accumulated jnean 

1 5 where accumulatedjnean was calculated in equation (8) and accumulated NLR Wiis calculated 
in equation (9). This value represents a weighted average of the net leg revenues 142 for each 
demand point from the top of the demand point list 140 up through the current demand point. 
A protection level 1 56 is determined for the accumulated demand point using a function (block 
133) further described hereinbelow in FIG. 1 5. Finally, the value of the next net leg revenue 1 57 

20 is set to the value of the net leg revenue 1 42 of the next demand point in the demand point list 
140. If there are no more demand points in the demand point list 140, the value of the next net 
leg revenue 1 57 is set to 0. If all of the demand points in the demand point list 140 have been 
processed (block 135), the marginal value 34 is set to the net leg revenue 142 of the last demand 
point in the demand point list 140 (block 136) and the function returns. Otherwise, if more 

25 demand points remain (block 135) and the difference of supply 44 minus protection level 155 
is greater than 0 (block 137), control continues at the top of the control loop (blocks 132-137) 
for processing of the next demand point. 

Otherwise, if the difference of supply 44 minus protection level 1 56 is less tlian or equal 
to 0 (block 137), the marginal value 34 is detemiined by using linear interpolation between the 

30 last two accumulated demand points (block 1 38) as expressed by the equation: 

(PL - - supply) {supply - PL - ) 

.ff . X [ + NNLRi X LL = (12) 

(PLi - PL,._,) ' (PL^ - PL,..,) 

3 5 where NNLR^ and NNLR^^j are the next net leg revenues 1 57 of the last two accumulated demand 
points, supply is the supply of seats 44, PL, and PLy.y are the protection levels 1 56 for the last 
two accumulated demand points and MV is the marginal value 34. 
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The linear interpolation will be explained with reference to FIG. 13 which shows a graph 
of a flight leg departure protection level curve generated by the ftuiction of FIG. 12. The x-axis 
represents the next net leg revenue 157 expressed in dollars and the y-axis represents the 
protection level 1 56 expressed in seats. Accumulated demand points 161A, 161B, 161C, 161D, 
161E are graphed as a function of their respective next net leg revenue 157 points 160A, 160B, 
160C, 160D, 160E. 

By way of example, assume that a marginal value 34 for the ORD-JFK flight leg departure 
31 is being determined with a supply of seats 44 for this flight leg departure equalling 10. 
Protection levels are calculated as illustrated, by way of example, in the table of FIG. 14C. Each 
row 151 A, 15 IB, 151C, 15 ID corresponds to an accumulated demand point as determined 
during each successive iteration of the control loop (blocks 132-137). Each coliunn corresponds 
to an accumulated mean 152, accumulated net leg revenue 153, accumulated variance 154, 
average net leg revenue 155, protection level 156 and next net leg revenue 157. Row 151 A 
represents the initialized value of the accumulated demand point. Row 151 B represents values 
determined during the first iteration of the control loop (block 132-137), row 151C during the 
second iteration and row 1 5 1 D during the third iteration. Using row 1 5 1 C as an example, the 
accumulated mean 152 is calculated using equation (8) as follows: 

2+4=6 

The accumulated net leg revenue 1 53 is calculated using equation (9) as follows: 

(365 X 2) + (322 x 4) = 2018 
The accumulated variance 1 54 is calculated using equation (1 0) as follows: 

1+3=4 

The average net leg revenue 155 is cdculated using equation (1 1) as follows: 

Mi . 336 

6 

The protection level 156 is determined by a function (block 134) which is described further 
hereinbelow in FIG. 1 5. The next net leg revenue 1 57 is 265 which is the value of the net leg 
revenue 142 of the third demand point 141 C on the demand point list 40, In this example, the 
result of the protection level calculation is 7. 

Starting from the first accumulated demand point 151 A, protection levels 156 are 
calculated and compared to supply 44. The first accumulated demand point 151 A has a 
protection level 1 56 of 0 which is less than the supply 44 of 10. This first point is plotted as 
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point 161 A in the graph of FIG. 13. The second accumulated demand point 151 A has a 
protection level 156 of 3 which is also less than the supply 44 of 10. This second point is plotted 

5 as point 1 6 IB in the graph of FIG. 13. The third accumulated demand point 151C has a 
protection level 156 of 7 v^ich is also less than the supply 44 of 10. This third point is plotted 
in the graph of FIG. 13 as point 161C. The fourth accumulated demand point 15 ID has a 
protection level 156 of 13 which is larger than the supply 44 of 10 which is the first demand 
point at which the difference of supply 44 minus protection level 1 56 is less than 0 (block 137). 

10 The difference of supply 44 minus protection level 156 is 10 - 1 3 = -3. Since -3 is less than 0 
(block 137), an interpolated marginal value 34 is determined between the accumulated demand 
points 161C and 161D by equation (12) as follows: 

$265 X (ILZJ^ . $205 X 00-7) ^ 
J5 (13 - 7) (13 - 7) 

Thus, the marginal value 33 is $235. The intersection 162 between the supply 44 (as indicated 
by line 1 63) and the interpolated demand (as indicated by line 164) falls at the point 165 on the 
X-axis where NNLR equals 235, the interpolated marginal value 34. Once this marginal value 
has been found (using equation (12)), the function retimis. 

20 FIG. 14A is a flow diagram of a demand point list creation function (block 130) as used 

by the function of FIG. 12. A flight path departure 36 containing the current flight leg departure 
31 is found by using the pointers 45 torn the flight path departures 36 to the flight leg departures 
3 1 (shown in FIG. 3). The marginal values 34 for all of the flight leg departures 3 1 contained 
in the flight path departure 36 except the cun^ent flight leg departure 3 1 are accumulated into a 

25 temporary local variable (block 1 5 1 ) as expressed by the equation: 

T^j (13) 

where MVj is the marginal value 34 for a flight leg departure 31 contained in the flight path 

30 departure 36 except the current flight leg departure 3 1 . 

When a marginal value for a particular flight leg departure 3 1 is being determined, the 
EMSR optimization function takes the net revenue, that is, booking class price 40, £md subtracts 
off the marginal values for all other flight leg departures 3 1 occurring in the current flight path 
departing 36 to obtain the net leg revenue 142 for the particular flight leg departure 31 being 

35 optimized. The net leg revenue 142 for the next point on the demand curve for the current flight 
leg departure 3 1 is determined (block 1 52) as expressed by the equation: 

BCPi-j:MVj (14) 
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where BCP^ is the booking class price 40 for the current flight path departure 36 and MVj is the 
marginal value 34 for a flight leg departure 31 contained in the flight path departure 36 except 
the current flight leg departure 31 . The net leg revenue 142 is stored in the demand point list 140 
as further described hereinbelow. 

If the net leg revenue 142 is greater than 0 (block 1 53), a new demand point is created and 
inserted into the demand point list 140 which is sorted in decreasing net leg revenue order (block 
1 54). For example, demand points for booking class prices 40 of $500, 400 and $ 1 50 are shown 
in FIG. 14B as having been inserted into the demand point list 140. Otherwise, if the net leg 
revenue is less than 0 (block 153), the current iteration over the points on the demand curve is 
terminated since any further points would also have net leg revenues less than 0. For example, 
demand points for booking class prices 40 of $ 1 30 and $ 1 20 will not be inserted into the demand 
point list 140. 

If the net leg revenue is still greater than 0 and more demand points corresponding to 
booking class prices 40 remain (block 155), control returns to the top of the inner control loop 
(block 1 52-1 55) for processing of the next demand point. Otherwise, if no further demand points 
remain (block 155) and more flight path departures 36 remain (block 156), control returns to the 
top of the outer control loop (blocks 1 50- 1 56) for processing of the next flight path departure 36. 
This process will continue to add points to the demand point list 1 40. If all flight path departures 
have been processed (block 1 56), the function returns. 

no. MB is a data strucmrc for a demand point list 140 as generated by the demand point 
list creation function of FIG. 14A. One demand point list 140 is created for each flight leg 
departure 3 1 that is to be optimized. The data structure for the demand point list 140 comprises 
demand point records 141 A, 141 B, 141 C, 141D, 141 E, 141 F containing net leg revenue 142, 
mean 143 and variance 144 fields. The net leg revenue 142 field stores the net leg revenue for 
the demand point corresponding to a booking class price 40 as calculated with equation (14) 
(block 152). The mean 41 and variance 42 for each demand record 39 are copied into the 
danand point list 140 as mean 143 and variance 144, respectively. For convenience, the net leg 
revenue 142 field will be referenced whenever a particular net leg revenue is referred to 
hereinbelow. 

FIG. 14B shows, by way of example, a demand point list 140 for the ORD-JFK flight leg 
departure 31. The first entry in this list, entry 141 A, corresponds to the first demand record 39 
in the demand records list 38 for the flight path departure 10 ORD-JFK-LHR-CDG as shown in 
FIG. 3. The first demand point encountered in the associated set of demand records 39 
corresponds to the booking class price 40 of $500. Using equation (14), the net leg revenue 142 
is: 
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BCPi - MVjjTK.^, - ^LHR-cDG = 



$500 - $70 - $65 = $365 



where £CP j is the first booking class price 40 of $500 for the ORD-JFK-LHR-CDG flight path 
departure 36 and MVjpjr^jjjj^ and MV^n_^[^ are the marginal values 34 for the J FK-LHR and 
LHR-CDG flight leg departures 31, respectively. Entries 141 C and 14IF are also Uiken from the 
flight path departure 10 ORD-JFK-LHR-CDG as shown in FIG. 14B. The other entries in the 
demand point list 140, entries 141B, 141D, 141E, are taken fix)m other flight path departures 
containing the ORD-JFK flight leg departure 10. The net leg revenue 142 value of $365 is stored 
in the demand point record 141 A. Similar calculations apply to the remaining demand point 
records 141B, 141C, 141D, 141E. 

FIG. 1 5 is a flow diagram of a protection level determination function (block 1 34) as used 
by the function of FIG. 12. Its purpose is to determine a protection level 156 for the current 
flight leg departure 31 being optimized as a function of the accumulated demanii point. The 
square root of the accumulated variance 154 is taken and stored in a local variable stddev (block 
1 80). The cumulative inverse normal of: 



is determined and stored in a local variable AT (block 181). ANLR is the average net leg revenue 
155 and NNLR is the next net leg revenue 157 of the accumulated demand point. In the 
described embodiment, the cumulative inverse normal is calculated using a numerical 
approximation, such as described in C. Hastings, Approximations for Digital Computers, p. 1 92, 
Princeton University Press (1955), the disclosure of which is incorporated herein by reference. 
The protection level 156 is determined (block 1 82) as expressed by the equation: 



\^4lere M is the accumulated mean 1 52 of the accumulated demand point and N and stdjdev are 
defined as described above. The function then returns. 

rv. Generic Mat^inal Value Determination Method Using Class Libraries 

As a further embodiment, the MVS 16 can be stmctured as a generic class library using 
an object-oriented programming language. In the described embodiment, the C-h- programming 
language is used; however, any object-oriented programming language can be employed. The 



(1 - 



NNLR 
ANLR 



PL ^ M ^ N ^ std dev 



(15) 
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discussion herein adopts the terminology used by the C-h- programming language, but is not 
meant to be restrictive as to the underlying concepts. 

FIG. 16 is a class diagram of a class library 189 for a method for determining marginal 
values for use in a perishable resource revenue management system. The class library 189 is 
structured into 3 layers: a generic layer 190, an industry-specific layer 191 and a client-specific 
layer 192. The generic layer 190 is used for determining marginal values for generic resources. 
In the described embodiment, a perishable resource problem is solved. The industry-specific 
layer 191 creates a set of derived classes for a particular industiy. In the described embodiment, 
derived classes for an airline, a hotel and a rental car agency are used, although other industries 
are feasible. The client-specific layer 192 creates a set of derived classes for a particular client 
relating to the corresponding industry-specific class. 

The generic layer 190 defines a base abstract resource class 1 93 and a perishable resource 
management problem (RMP) class 1 94. The resource class defines an interface for derived 
classes in the industry-specific layer 191 for a flight leg departure resource class 195 A, a room 
night resource class 195B and a rental car day resource class 195C. Other derived resource 
classes are feasible. The resource class 193 comprises two pure virtual fimctions supplyQ and 
demandQ (not shown). As further described hereinbelow, the ftinctions defined in portions of 
the class definitions for the industry-specific layer 191 and the client-specific layer 192 override 
the pure virtual fimctions defined in the generic layer 190. 

The perishable RMP class 194 is also in the generic layer 190 and defines an interface 
for derived abstract classes in the industry-specific layer 191 for solving an airline resource 
management problem 196 A, a hotel resource management problem 196B and a rental car 
resource management problem 196C. Other derived abstract classes are feasible. The line 
connecting the perishable RMP 1 94 with the resource class 1 93 indicates that the RMP has many 
resources to manage. 

The client-specific layer 192 comprises a set of derived classes defining the particular 
data manipulation fimctions used by a particular client's implementation to load and save data 
in the client's revenue management system 15 and to map the data into the generic data 
stmctuies. For instance, the airline #1 resource management problem class 197 A contains those 
loading, saving and mapping fimctions needed to utilize the generic data structures and methods 
for determining marginal values for airline #1 . This likewise applies to the airline #2 resource 
management problem 197B, the hotel #1 resource management problem 198 A, the hotel Wl 
resource management problem 198B, the company #1 resource management problem 199A and 
the company #2 resource management problem 199B classes. Other derived classes under 
industry-specific layer 191 class are feasible. 



-25- 



wo 97/25< 




PCT/US97/00738 



The resource class 1 93 defines a base abstract class for a generic resource for which 
marginal values will be detennined. The particular type of resource is defined in the derived 

5 classes defined in the industry-specific layer 191. There is a one-to-one correlation between the 
classes derived bom the resource class 193 and the specific abstract classes derived fix)m the 
perishable RMP class 194. The flight leg resource class 1 95 A corresponds to the airline resource 
management problem class 196A. The same applies to the room night resource class 195B and 
the hotel resource management problem class 196B and the rental car resource class 195C and 

10 the rental car resource management problem class 196C. The particular type of resource being 
solved is germane to the industry-specific layer 191 and consequently, it is unnecessary to define 
each resource abstractly. In other words, a flight leg departure is, for example, the common 
resource used by the airline resource management problem 196A irrespective of the particular 
airline. 

1 5 FIG. 1 7 is a class diagram of an instantiation 200 of the class library 1 89 of FIG. 1 6 for 

an airline reservations system. The base abstract perishable RMP class 194 defines a generic 
fiinction solveQ which uses the supplyQ and demandQ fimctions of the resource class 1 93 to 
solve for generic marginal values. This abstract class has a derived airline RMP class 196A 
which defines a set of data structures as shown in FIG. 3 for an airline-specific implementation. 

20 The airline RMP class 1 96 A has a derived airline #1 RMS class 1 97A which contains two airline 
#1 specific functions, loadQ and storeQ. These fiinctions embody the necessary functionality 
for loading flight leg, flight path and passenger value data from and storing marginal values to 
airline #rs airline revenue management system 15 (shovm in FIG. 2 A). 

The solveQ function of the RMP class 194 is further described hereinbelcw beginning 

25 in FIG. 18, The resource class 193 contains a pair of pure virtual functions supplyQ and 
demandQ that are overridden by the particular supplyQ and demandQ routines used in the 
instantiated problem, that is, by the funcrions contained in the particular industry-specific layer 
191 class. 

FIG. 1 8 is a flow diagram of a method for generically determining marginal values for 
30 use in a client-specific perishable resource management system. Data is loaded using the client- 
specific loadQ function (block 210). As described hereinabove, the loadQ function is defined 
in the client-specific layer 192 as a derived class member for a particular client. The loaded data 
is mapped into the problem space for the industry-specific layer 191 and generic layer 1 90 (block 
211). This mapping is accomplished through inheritance whereby the parent and grandparent 
35 classes, that is, the industry-specific layer 191 derived abstract classes and generic layer 190 
abstract class, instantiate client-specific objects to industry-specific and finally generic objects. 
Generic marginal values are then determined using the generic solveQ function, such as shown 
in the perishable RMP class 194 (block 212) as further described hereinbelow in FIG, 19. The 



-26- 



wo 97/2! 



PCT/US97/00738 



generic marginal values are mapped back into the client-specific problem space, that is, to the 
particular client-specific layer 192 class (block 213). Finally, the marginal values as mapped 
fi-om generic marginal values to client-specific marginal values are stored using the storeQ 
fimction contained in the client-specific layer 192 class (block 214). 

FIGS. 19-26 are flow diagrams of an iterative fimction for determining generic marginal 
values for use in a perishable resource management system. The opemtions described in these 
figures are the same as those previously described for an airline-specific implementation, 
although the two implementations operate in different problem spaces. In particular, the program 
variables used by the generic class library 1 89 are presented in terms of perishable resources 
rather than seats on flight leg departures. The correlation between the airline-specific fimctions 
and the generic class library fimctions are shown below in Table II. FIGS. 10 and 23 are 
identical, so there is no entry in Table II for these figures. Note that identical reference numerals 
are used for the flow diagrams in FIGS. 5-6, 8-10, 12 and 14A-15 as in FIGS. 19-26, 
respectively, except that an apostrophe sign has been added to the reference numerals in the 
latter. Therefore, it will be understood that the same portions of the Detailed Description with 
respect to FIGS. 5-6, 8-10, 12 and 14A-15 (pertaining to the airiine-specific method) also apply 
to FIGS. 19-26 (pertaining to the generic class library method) with exceptions as noted in Table 
n for program variable names. 

Table II 
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Table II 
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While various embodiments of the invention have been particularly shown, it will be 
understood by those skilled in the art that the various changes in form and detail may be made 
without departing bom the scope and spirit of the present invention, wherein: 

30 
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WHAT IS CLAIMED IS: 

1 . A method using a marginal value system (16) for determining marginal values 
(34) for perishable resources (3 1 ), comprising the steps of: 

loading data (50) for the perishable resources (3 1 ) and composite resources (36) 
fiom a perishable resource revenue management system (15) into the marginal value system ( 1 6); 

constructing (51) internal data structures (30, 35, 38) for linking each of the 
perishable resources (3 1) to their associated composite resources (36) and for linking each of the 
composite resources (36) to their associated perishable resources (31); 

detemiining (52) the marginal values (34) for the perishable resources (3 1 ) using 
a continuous optimization function (FIGS. 20 and 24) using interdependencies among the 
perishable resources (31) and the composite resources (36) in the internal data structures (30, 35, 
38); and 

storing (53) the maiginal values (34) from the marginal value system (16) into the 
perishable resource revenue management system (1 5). 

2. A method according to claim 1, the step of determining further comprising the 
steps of: 

evaluating (62) a locally optimal marginal value (34) for the perishable resources 
(31) using the continuous optimization function (FIGS. 20 and 24) on each of the perishable 
resources (31); and 

performing (66) the step of evaluating (62) in successive iterations until a globally 
optimal marginal value (34) is evaluated for the perishable resources (31) based on a 
convergence criterion (67). 

3. A method according to claim 2, further comprising the steps of: 
monitoring (63, 64, 65) a maximum difference between such locally optimal 

marginal values (34) evaluated during a current iteration and such locally optimal marginal 
values (34) evaluated during a previous iteration; 

comparing (67) the maximum difference and the convergence criterion upon 
completion of the current iteration whereby convergence is satisfied when the maximum 
difference is less than the convergence criterion. 
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4. A method according to claim I, the step of determining further comprising the 
steps of: 

evaluating (62) a locally optimal marginal value (34) for one of the perishable 
resources (31) using the continuous optimization function (FIGS. 20 and 24); and 

perfomfiing (66) the step of evaluating (62) in successive iterations for another of 
the perishable resources (31) until a locally optimal marginal value (34) has been evaluated for 
all of the perishable resources (3 1 ). 

5. A method according to claim 4, further comprising the steps of: 
determining (62) the absolute difference between such locally optimal marginal 

value (34) evaluated for one of the perishable resources (31) during a current iteration and such 
locally optimal marginal value (34) evaluated for the resource during a previous iteration; 

comparing (63) the absolute difference to a maximum difference between such 
locally optimal marginal value (34) evaluated for another of the perishable resources (31) during 
the current iteration and such locally optimal marginal value (34) evaluated for the resource 
during the previous iteration whereby the maximum difference equals or exceeds all such 
absolute differences between such locally optimal marginal values (34) evaluated for all other 
perishable resources (31) during a current iteration and such locally optimal marginal values (34) 
evaluated for the resources during a previous iteration; and 

replacing (65) the maximum difference with the absolute difference when the 
absolute difference is greater than the maximiun difference. 

6. A method according to claim 1 , further comprising the steps of: 

loading (50) resource value data comprising a supply (44) of the perishable 
resources (31) and a set of resource prices (40) from the perishable resource revenue 
management system (15) into the marginal value system (16); and 

constmcting (51) an internal data stmcture (38) for linking each of the composite 
resources (36) to their associated set of resource prices (40) and for linking each of the perishable 
resources (31 ) to their associated supply (44). 

7. A method according to claim 6, whereby the continuous optimization function is 
a deterministic optimization function, the step of determining (52) further comprising evaluating 
(62) locally optimal marginal values (34) based on a certain demand (81, 85) for each of the 
composite resources (36). 
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8. A method according to claim 7, whereby the detenninistic optimization function 
is a supply-demand balance optimization fimction (FIG. 20), fiirther comprising the step of: 

searching (73) for a marginal value (34) for one of the composite resources (36) 
whereby the supply (44) substantially equals the demand for the one composite resource (36). 

9. A method according to claim 8, further comprising the steps of: 
determining (71) the difference between the supply (44) and the certain demand 

(81, 85) for the one resource (36); and 

setting (70) the marginal value (31) for the one resource (36) to indicate that the 
supply (44) exceeds the demand when the difference is positive (72). 

1 0. A method according to claim 8, further comprising the step of: 
detemiining (71, 73) the certain demand (81, 85) for the one perishable resource 

(31) based on the marginal values (34) for each other of the perishable resources (3 I ) contained 
in each of the composite resources (36) containing the one perishable resource (31 ). 

11. A method according to claim 1 0, further comprising the steps of: 
determining (91) a composite resource demand for each of the composite 

resources (36) containing the one perishable resource (31); and 

adding (92) the composite resource demand for each of the composite resources 
(3 6) to the certain demand (8 1 , 85) for the one perishable resource (31). 

12. A method according to claim 1 1 , further comprising the steps of: 

adding (101) the marginal value (34) for each of the perishable resources (31) 
contained in one of the composite resources (36) to a total marginal value; and 

detennining (103) the composite resource demand for the one composite resource 
(36) using the total marginal value. 

13. A method according to claim 1 2, whereby each of the set of resource prices (40) 
comprises a comesponding demand point (120 A, 120B, 120C, 120D, 120E) on a demand curve, 
further comprising the steps of: 

searching (110) for the first demand point on the demand curve (120 A, 120B, 
120C, 120D, 120E) conesponding to a resource price (40) that exceeds the total marjpnal value; 

detennining (1 12) a linear average (123) between the first demand point (120D) 
and a second demand point (I20C) previous to the first demand point (120D); and 



-32- 



wo 97/25( 



PCT/US97/00738 



deterrniniiig (1 1 3) the composite resource demand based on an intersection of the 
supply (44) and the linear average (123). 

5 

1 4. A method according to claim 6, whereby the continuous optimization function is 
a non-deterministic optimization function, the step of determining (52) further comprising 
evaluating (62) locally optimal maiginal values (34) based on an uncertain demand ( 1 6 1 A, 1 6 1 B, 
161C, 16 ID) for each of the composite resources (36). 

10 

15, A method according to claim 1 4, whereby the continuous optimization function 
is a nonnletenninistic optimization function based on a supply of each of the perishable resources 
(31) and an uncertain demand for the composite resources. 

15 16. A method according to claim 15, whereby the non-deterministic optimization 

function is an expected marginal resource revenue (EMRR) optimization function (FIG. 24), 
further comprising the steps of; 

searching (138) for a mai^inal value (34) for one of the composite resources (36) 
whereby the supply (44) substantially matches a protection level (162) for the one composite 

20 resource (36). 

1 7. A method according to claim 1 6, fiirther comprising the step of: 

loading (50) resource value data further comprising a set of means (4 1 ) and a set 
of variances (42) from the perishable resource revenue management system (15) into the 
25 marginal value system (16) whereby the set of means (41) and the set of variances (42) 
correspond to the set of resource prices (40); and 

constructing (51) an internal data structure (38) for hnking each of the composite 
resources (36) to their associated set of means (41) and set of variances (42). 

30 18. A method according to claim 1 6, further comprising the steps of: 

creating (130) a demand point list (140) for the one perishable resource (31) 
containing a set of demand points (160A, 160B, 160C. 160D) corresponding to each of the set 
of resource prices (40); 

determining (132) an accumulated mean (152), an accumulated net leg revenue 
35 (1 53) and an accumulated variance based on the set of resource prices (40), the set of means (4 1 ) 
and the set of variances (42) for each demand point in the demand point list (140); 

determining (133) an updated average net resource revenue (154) based on the 
accumulated net leg revenue (1 53) and the accumulated mean (1 52); and 
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determining (134) the protection level (162) for each of a set of accumulated 
demand points ( 1 6 1 A, 1 6 1 B, 1 6 1 C, 1 6 1 D) using the accumulated mean ( 1 52), the accumulated 
variance and the average net leg revenue (154) for each demand point (160A, 160B, I60C, 
160D) in the demand point list (140). 

19. A method according to claim 18, whereby the demand point list (140) is 
determined, fiuther comprising the steps of: 

accumulating (151) the marginal values (34) for all of the perishable resources 
(31) in a composite resource (36) containing the one perishable resoim:e (31) except for the 
marginal value (31) for the one composite resource (34) for each of the set of resource prices 
(40); 

determining (1 51 ) a net leg revenue (142) for each demand point (160A, 1 603, 
160C, 160D) in the demand point list (140) using the accumulated marginal values (34); and 

creating (154) a new demand point (160 A, 160B, 160C, 160D) when the net leg 
revenue (142) is positive (153). 

20. A method according to claim 1 9, further comprising the steps of: 

setting (152) the net leg revenue (142) to the difference of the corresponding 
resource price (40) minus the accumulated marginal values. 

21. A method according to claim 1 8, further comprising the steps of: 
determining (132) the accumulated mean (152) comprising a summation of each 

mean in the set of means (4 1 ) for each of the set of demand points ( 1 60A, 1 608, 1 60C, 1 60D); 

determining (132) the accumulated net leg revenue (1 53) comprising a summation 
of each of the set of resource prices (40) multiplied by each of the means in the set of means (4 1 ) 
for each of the set of demand points (160A, 160B, 160C, 160D); and 

determining (132) the accimiulated variance comprising a summation of each 
variance in the set of variances (41) for each of the set of demand points (160 A, 160B, 160C, 
160D). 

22. A method according to claim 1 8, further comprising the steps of: 
determining (133) the updated average net resource comprising the quotient of 

the accumulated net leg revenue ( 1 53) divided by the accumulated mean (1 52). 
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23 . A method according to claim 1 8, further comprising the steps of: 
determining (1 81) a standard deviation equalling a square root of the accumulated 

variance; 

determining (183) an inverse cumulative normal of the equation: 
(1-NNRR/ANRR) 

whereby NNRR corresponds to next net resource revenue ( 1 57) for a demand point ( 1 60 A, 160B, 
160C, 160D) in the demand point list (140) and ANRR corresponds to the average net resource 
revenue (154) for each demand point (160A, 160B, 160C, 160D) in the demand point list (140); 
and 

setting (1 85) the protection level ( 1 62) to the sum of the accumulated mean ( 1 52) 
plus the inverse cumulative normal multiplied by the standard deviation, 

24. A marginal value system (16) for determining marginal values (34) for perishable 
resources (3 1 ), comprising: 

means for loading data (50) for the perishable resources (3 1 ) and composite 
resources (36) from a perishable resource revenue management system (15) into the marginal 
value system (16); 

means for constructing (51) internal data structures (30, 35, 38) for linking each 
of the perishable resources (3 1) to their associated composite resources (36) and for linking each 
of the composite resources (36) to their associated perishable resources (31); 

means for detennining (52) the marginal values (34) for the perishable resoxuces 
(31) using a continuous optimization function (FIGS. 20 and 24) using the perishable resources 
(3 1 ) and the composite resources (36) in the internal data structures (30, 35, 38); and 

means for storing (53) the marginal values (34) from the marginal value system 
(16) into the perishable resource revenue management system (15). 

25. A method of producing a marginal value (34) representing currency determined 
using a marginal value system (16) for use in a perishable resource revenue management system 
(14) for granting and denying (22, 23) a sale of one or more perishable resources (31) expiring 
at a future time (45) depending on marginal values (34) for each of such perishable resources 
(31) received from a perishable resource revenue management system (36), comprising the steps 
of: 

loading data (50) for the perishable resources (31), composite resources (36) and 
a resource value (40, 41, 42) from a perishable resource revenue management system (15) into 
a marginal value system (16); 
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constructing (51) internal data structures (30, 35, 38) for the perishable resources 
(31), composite resources (36) and resource value (40, 41 , 42); 

determining (52) marginal values (34) for the perishable resources (3 1 ) using a 
continuous optimization function (FIGS. 20 and 24); and 

storing (53) the marginal values (34) fix)m the marginal value system ( 1 6) into the 
perishable resource revenue management system (15). 

26. A method using a computer (15) for determining marginal values (34) for 
perishable resources (3 1 ), comprising the steps of: 

evaluating (62) a locally optimal marginal value (34) for one of the perishable 
resources (31) using a continuous optimization function (FIGS. 20 and 24) dependent on the 
marginal values (34) for the other perishable resources (3 1); and 

iteratively reevaluating (67) the locally optimal marginal value (34) until a 
globally optimal marginal value (34) is attained. 

27. A method according to claim 26, further comprising the steps of: 
iteratively performing (66) the step of evaluating (62) until a locally optimal 

marginal value (34) is attained for each of the perishable resources (31). 

28. A method according to claim 26, whereby the continuous optimizaiion function 
(FIGS. 20 and 24) is a supply demand balance optimization function (FIG. 20), further 
comprising the steps of: 

determining (72) whether a supply (80, 81) for the one perishable resource (31) 
minus a demand (81, 85) for the one perishable resource (31) is positive; and 

searching (73) for the locally optimal marginal value (34) that makes the supply 
(80, 8 1 ) equal the demand (81,85). 

29. A method according to claim 26, whereby the continuous optimization function 
(FIGS. 20 and 24) is a expected marginal resource revenue optimization function (FIG. 24), 
further comprising the steps of: 

creating (1 30) a demand point list (140) for the one perishable resource (31) based 
on a set of resource prices (40), means (41) and variances (42) associated with each of the 
composite resources (36) containing the one perishable resource (31); 

determining (134) a protection level (162) for each demand point (1 60 A, 160B, 
160C, 160D) in the demand point list (140); and 
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searching (138) for the locally optimal marginal value (34) that makes the supply 
(80, 81) equal the protection level (162). 

30. A marginal value system (16) for determining marginal values (34) for perishable 
resources (31), comprising: 

means for evaluating (62) a locally optimal marginal value (34) for one of the 
perishable resources (31) using a continuous optimization fimction (FIGS. 20 and 24) dependent 
on the marginal values (34) for the other perishable resources (3 1 ); and 

means for iteratively reevaluating (67) the locally optimal marginal value (34) 
until a globally optimal marginal value (34) is attained. 

31. A method using a marginal value system (16) for determining marginal values 
(34) for seats on flight leg depamires (31) departing on a future departure date (45), comprising 
the steps of: 

loading data (50) for flight leg departures (31) and flight path departures (36) 
from an airline revenue management system (15) into the marginal value system (16); 

constructing (51) internal data structures (30, 35, 38) for linking each of the flight 
leg departures (3 1) to their associated flight path departures (36) and for linking each of the flight 
path departures (36) to their associated flight leg departures (31); 

determining (52) the marginal values (34) for seats on the flight leg departures 
(31) using a continuous optimization function (FIGS. 5 and 12) using interdependencies among 
the flight leg departures (31) and the flight path departures (36) in the internal data structures 
(30, 35, 38); and 

storing (53) the marginal values (34) from the marginal value system (1 6) into the 
airline revenue management system (15). 

32. A method according to claim 3 1 , the step of determining further comprising the 
steps of: 

evaluating (62) a locally optimal marginal value (34) for seats on the flight leg 
departures (31) using the continuous optimization function (FIGS. 5 and 12) on each of the flight 
leg departures (3 1 ); and 

performing (66) the step of evaluating (62) in successive iterations until a globally 
optimal marginal value (34) is evaluated for the flight leg departures (31) based on a 
convergence criterion (67). 
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33. A method according to claim 32, further comprising the steps of: 
monitoring (63, 64, 65) a maximum difference between such locally optimal 

marginal values (34) evaluated during a current iteration and such locally optimal marginal 
values (34) evaluated during a previous iteration; 

comparing (67) the maximum difference and the convergence criterion upon 
completion of the current iteration whereby convergence is satisfied when the maximum 
difference is less than the convergence criterion. 

34. A method according to claim 31, the step of detennining further comprising the 
steps of: 

evaluating (62) a locally optimal marginal value (34) for a seat on one of the flight 
leg departures (3 1 ) using the continuous optimization function (FIGS. 5 and 1 2); and 

performing (66) the step of evaluating (62) in successive iterations for another of 
the flight leg departures (31) until a locally optimal marginal value (34) has been evaluated for 
all of the flight leg departures (31). 

35. A method according to claim 34, further comprising the steps of: 
determining (62) the absolute difference between such locally optimal marginal 

value (34) evaluated for a seat on one of the flight leg departures (31) during a current iteration 
and such locally optimal marginal value (34) evaluated for the seat during a previous iteration; 

comparing (63) the absolute difference to a maximimi difference between such 
locally optimal marginal value (34) evaluated for a seat on another of the flight leg departures 
(31) during the current iteration and such locally optimal marginal value (34) evaluated for the 
seat during the previous iteration whereby the maximiun difference equals or exceeds all such 
absolute differences between such locally optimal marginal values (34) evaluated for seats on 
all other flight leg departures (31) during a current iteration and such locally optimal marginal 
values (34) evaluated for the seats during a previous iteration; and 

replacing (65) the maximimi difference with the absolute difference when the 
absolute difference is greater than the maximum difference. 

36. A method according to claim 3 1 , further comprising the steps of: 

loading (50) passenger value data comprising a supply of seats (44) on each of 
the flight leg departures (31) and a set of booking class prices (40) from the airline revenue 
management system (15) into the marginal value system (16); and 
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constructing (51) an internal data structure (38) for linking each of the flight path 
departures (36) to their associated set of booking class prices (40) and for linking each of the 
5 flight leg departures (3 1 ) to their associated supply of seats (44). 

37. A method according to claim 36, whereby the continuous optimization ftinction 
is a deterministic optimization function, the step of determining (52) further comprising 
evaluating (62) locally optimal marginal values (34) based on a certain demand (8 1 , 85) for seats 

1 0 on each of the flight leg departures (36). 

38. A method according to claim 37, whereby the deterministic optimization function 
is a supply-demand balance optimization function (FIG. 5), further comprising the step of: 

searching (73) for a marginal value (34) for one of the flight leg departures (36) 
1 5 whereby die supply of seats (44) substantially matches the demand for seats on the one flight leg 
departure (36). 

39. A method according to claim 38, further comprising the steps of 
determining (71) the difference between the supply of seats (44) and the certain 

20 demand (8 1 , 85) for seats on the one flight leg departure (36); and 

setting (70) the marginal value (31) for the one flight leg departure (36) to 
indicate that the supply of seats (44) exceeds the demand for seats when the diflference is positive 
(72). 

25 40. A method according to claim 38, further comprising the step of: 

determining (71, 73) the certain demand (81 , 85) for the one flight leg departure 
(3 1 ) based on the marginal values (34) for each other of the flight leg departures (3 1 ) contained 
in each of the flight path departures (36) containing the one flight leg departure (3 1 ). 

30 41 . A method according to claim 40, further comprising the steps of: 

determining (91) a path demand for each of the flight path departures (36) 

containing the one flight leg departure (3 1 ); and 

adding (92) the path demand for each of the flight path departures (36) to the 

certain demand (81 , 85) for the one flight leg departure (3 1). 

35 

42. A method according to claim 41 , further comprising the steps of: 

adding (101) the marginal value (34) for each of the flight leg departures (3 1) 
contained in one of the flight path departures (36) to a total marginal value; and 
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determining (103) the path demand for the one flight path departure (36) using 
the total marginal value. 

5 

43. A method according to claim 42, whereby each of the set of booking class prices 
(40) comprises a corresponding demand point (120 A, 120B, I20C, 120D, 120E) on a demand 
curve, further comprising the steps of: 

searching (1 10) for the first demand point on the demand curve (120A, 120B, 
10 I20C, 1 20D, 120E) corresponding to a booking class price (40) that exceeds the total marginal 
value; 

determining (1 12) a linear average (123) between the first demand point (1 20D) 
and a second demand point (120C) previous to the first demand point (120D); and 

deteimining (1 1 3) the path demand based on an intersection of the supply of seats 
1 5 (44) and the linear average (1 23). 

44. A method according to claim 36, whereby the continuous optimization function 
is a non-deterministic optimization function, the step of determining (52) further comprising 
evaluating (62) locally optimal marginal values (34) based on an uncertain demand (161 A, 161 B, 

20 1 6 1 C, 1 6 1 D) for seats on each of the flight leg departures (36). 

45. A method according to claim 44, whereby the continuous optimization function 
is a non-deterministic optimization function based on a supply of seats (1 55) on each of the flight 
leg departures (3 1) and an uncertain demand for the seats. 

25 

46. A method according to claim 45, whereby the non-deterministic optimization 
function is an expected marginal seat revenue (EMSR) optimization function (FIG. 12), further 
comprising the steps of; 

searching (1 38) for a marginal value (34) for one of the flight leg departures (36) 
30 Avhereby the supply of seats (44) substantially matches a protection level (1 62) for seats on the 
one flight leg departure (36). 

47. A method according to claim 46, further comprising the step of: 

loading (50) passenger value data further comprising a set of means (41) and a 
35 set of variances (42) from the airline revenue management system (15) into the marginal value 
system (1 6) whereby the set of means (41) and the set of variances (42) correspond to the set of 
booking class prices (40); and 
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constructing (51) an internal data istructure (38) for linking each of the flight path 
departures (36) to their associated set of means (41) and set of variances (42). 

48. A method according to claim 46, fiirther comprising the steps of: 

creating (130) a demand point list (140) for the one flight leg departure (31) 
containing a set of demand points (160 A, 160B, 160C, 160D) corresponding to each of the set 
of booking class prices (40); 

determining (132) an accumulated mean (152), an accumulated net leg revenue 
(153) and an accumulated variance based on the set of booking class prices (40), the set of means 
(41) and the set of variances (42) for each demand point in the demand point hst (140); 

determining (133) an updated average net leg revenue (154) based on the 
accumulated net leg revenue (153) and the accumulated mean (152); and 

determining (134) the protection level (162) for each of a set of accumulated 
demand ix)ints( 161 A, 161B, 161C, I61D) using the accumulated mean (152), the accumulated 
variance and the average net leg revenue (154) for each demand point (160A, 160B, 160C, 
160D) in the demand point list (140). 

49. A method according to claim 48, whereby the demand point list (140) is 
determined, further comprising the steps of: 

accumulating (151) the margmal values (34) for all of the flight leg departures 
(3 1) in a flight path departure (36) containing the one flight leg departure (31) except for the 
marginal value (31) for the one flight leg departure (34) for each of the set of booking class 
prices (40); 

determining (151) a net leg revenue (142) for each demand point (160 A, 160B, 
1 60C, 1 60D) in the demand point list (140) using the accumulated marginal values (34); and 

creating ( 1 54) a new demand point ( 1 60 A, 1 60B, 1 60C, 1 60D) when the net leg 
revenue (142) is positive (153), 

50. A method according to claim 49, further comprising the steps of: 

setting (152) the net leg revenue (142) to the difference of the corresponding 
booking class price (40) minus the accumulated marginal values. 

51. A method according to claim 48, further comprising the steps of: 
determining (132) the accumulated mean (152) comprising a summation of each 

mean in the set of means (41)foreach ofthe set ofdemand points (160 A, 160B, 160C, 160D); 
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detenriining (1 32) the accumulated net leg revenue ( 1 53) comprising a summation 
of each of the set of booking class prices (40) multiplied by each of the means in the set of means 
5 (41) for each of the set of demand points (160A, 160B, 160C, 160D); and 

detemiining (132) the acciimulated variance comprising a summation of each 
variance in the set of variances (41) for each of the set of demand points (160A, 160B, 160C, 
160D). 

10 52. A method according to claim 48, fijrther comprising the steps of: 

delerminmg (133) the updated average net leg comprising the quotient of the 
accumulated net leg revenue (153) divided by the accumulated mean (152). 

53. A method according to claim 48, further comprising the steps of: 

1 5 determining (1 8 1 ) a standard deviation equalling a square root of the accumulated 

variance; 

determining (183) an inverse cumulative normal of the equation: 
(1 -NNLR/ANLR) 

whereby NNLR corresponds to a next net leg revenue (157) for a demand point (1 60 A, 160B, 
20 160C, 160D) in the demand point list (140) and ANLR corresponds to the average net leg 
revenue (154) for each demand point (160A, 160B, 160C, 160D) in the demand point list (140); 
and 

setting (185) the protection level (162) to the sum of the accumulated mean (152) 
plus the inverse cumulative normal multiplied by the standard deviation. 

25 

54. A marginal value system (16) for determining marginal values (34) for seats on 
flight leg departures (31) departing on a future departure date (45), comprising: 

means for loading data (50) for flight leg departures (31) and flight path 
departures (36) from an airline revenue management system (15) into the marginal value system 
30 (16); 

means for constmcting (51) intemal data structures (30, 35, 38) for linking each 
of the flight leg departures (3 1 ) to their associated flight path departures (36) and for linking each 
of the flight path departures (36) to their associated flight leg departures (3 1); 

means for determining (52) the marginal values (34) for seats on the flight leg 
35 departures (31) using a continuous optimization function (FIGS. 5 and 12) using the flight leg 
departures (3 1 ) and the flight path departures (36) in the intemal data structures (30, 35, 38); and 

means for storing (53) the marginal values (34) from the marginal value system 
(16) into the airline revenue management system (15). 
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55. A method of producing a marginal value (34) representing currency determined 
using a marginal value system (16) for use in conjunction with an airline reservations system 
(14) for granting and denying (22, 23) a sale of seats on one or more flight leg departures (31 ) 
departing on a future departure date (45) depending on marginal values (34) for seats on each 
of such flight leg departures (31) received from an airline revenue management system (36), the 
method comprising the steps of: 

loading data (50) for flight leg departures (31), flight path departures (36) and 
passenger value (40, 41, 42) from an airline revenue management system (15) into a marginal 
value system (16); 

constructing (51) internal data structures (30, 35,38) for the flight leg departures 
(3 1 X flight path departures (36) and passenger value (40, 41, 42); 

determining (52) marginal values (34) for seats on the flight leg departures (31) 
using a continuous optimization function (FIGS. 5 and 12); and 

storing (53) the marginal values (34) from the marginal value system (16) into the 
airline revenue management system (15). 

56. A method using a computer (1 5) for determining marginal values (34) for seats 
on flight leg depsirtures (31) departing on a future departure date (45), comprising the steps of: 

evaluating (62) a locally optimal marginal value (34) for a seat on one of the flight 
leg departures (3 1 ) using a continuous optimization function (FIGS. 5 and 1 2) dependent on the 
marginal values (34) for a seat on the other flight leg departures (31); and 

iteratively reevaluating (67) the locally optimal marginal value (34) until a 
globally optimal marginal value (34) is attained. 

57. A method according to claim 56, further comprising the steps of: 
iteratively performing (66) the step of evaluating (62) until a locally optimal 

margmal value (34) is attained for each of the flight leg departures (3 1 ). 

58. A method according to claim 56, whereby the continuous optimization function 
(FIGS. 5 and 12) is a supply demand balance optimization function (FIG. 5), further comprising 
the steps of: 

determining (72) whether a supply of seats (80, 8 1) for the one flight leg departure 
(3 1 ) minus a demand for seats (8 1 , 85) for the one flight leg departure (3 1 ) is positive; and 

searching (73) for the locally optimal marginal value (34) that makes the supply 
of seats (80, 81 ) equal the demand for seats (81, 85). 
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59. A method according to claim 56, whereby the continuous optimization function 
(FIGS- 5 and 12) is a expected marginal seat revenue optimization function (FIO. 12), further 
comprising the steps of; 

creating (1 30) a demand point list (140) for the one flight leg departure (3 1 ) based 
on a set of booking class prices (40), means (41) and variances (42) associated with each of the 
flight path departures (36) containing the one flight leg departure (31); 

determining (134) a protection level (162) for each demand point (1 60A, 160B, 
160C, 160D) in the demand point list (140); and 

searching (138) for the locally optimal marginal value (34) that makes the supply 
of seats (80, 81) equal the protection level (162). 

60. A marginal value system (16) for determining marginal values (34) for seats on 
flight leg departures (3 1 ) departing on a fiiture departure date (45), comprising: 

means for evaluating (62) a locally optimal marginal value (34) for a seat on one 
of the flight leg departures (31) using a continuous optimization function (FIGS. 5 and 12) 
dependent on the marginal values (34) for a seat on the other flight leg departures (31); and 

means for iteratively reevaluating (67) the locally optimal marginal value (34) 
until a globally optimal marginal value (34) is attained. 
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